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INTERACTIONS 

By D. Bertacchi, N. Lanchier^ and F. Zucca 
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We introduce spatially explicit stochastic processes to model mul- 
tispecies host-symbiont interactions. The host environment is static, 
modeled by the infinite percolation cluster of site percolation. Sym- 
bionts evolve on the infinite cluster through contact or voter type 
interactions, where each host may be infected by a colony of sym- 
bionts. In the presence of a single symbiont species, the condition for 
invasion as a function of the density of the habitat of hosts and the 
mciximal size of the colonies is investigated in details. In the pres- 
ence of multiple symbiont species, it is proved that the community of 
symbionts clusters in two dimensions whereas symbiont species may 
coexist in higher dimensions. 

1. Introduction. The term symbiosis was coined by the mycologist Hein- 
rich Anto de Bary to denote close and long-term physical and biochemical 
interactions between different species, in contrast with competition and pre- 
dation that imply only brief interactions. Symbiotic relationships involve 
a symbiont species, smaller in size, that always benefits from the relationship, 
and a host species, larger in size, that may either suffer, be relatively unaf- 
fected, or also benefit from the relationship, which are referred to as parasis- 
tism, commensalism, and mutualism, respectively. The degree of specificity 
of the symbiont is another important factor: while some symbionts may live 
in association with a wide range of host species, in which case the sym- 
biont is called a generalist, others are highly host-specific indicating that 
they can only benefit from few host species. Symbiotic relationships, either 
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pathogenic or mutualistic, are ubiquitous in nature. For instance, more than 
90% of terrestrial plants [28] live in association with mycorrhizal fungi, with 
the plant providing carbon to the fungus and the fungus providing nutrients 
to the plant, most herbivores have mutualistic gut fauna that help them 
digest plant matter, and almost all free-living animals are host to one or 
more parasite taxa [27]. 

To understand the role of spatial structure on the persistence of host- 
parasite and host-mutualist associations, Lanchier and Neuhauser [21-23] 
have initiated the study of multispecies host-symbiont systems including lo- 
cal interactions based on interacting particle systems. The stochastic process 
introduced in [22] describes the competition among specialist and generalist 
symbionts evolving in a deterministic static environment of hosts. The math- 
ematical analysis of this model showed that fine-grained habitats promote 
generalist strategies, while coarse-grained habitats increase the competitive- 
ness of specialists. The stochastic process introduced in [21, 23] includes in 
addition a feedback of the hosts, which is modeled by a dynamic-host system. 
This process has been further extended by Durrett and Lanchier [13]. The 
host population evolves, in the absence of symbionts, according to a biased 
voter model, while the symbiont population evolves in this dynamic envi- 
ronment of hosts according to a contact type process. The parameters of the 
process allow to model the effect of the symbionts on their host as well as the 
degree of specificity of the symbionts, thus resulting in a system of coupled 
interacting particle systems, each describing the evolution of a trophic level. 
The model is designed for the understanding of the role of the symbionts 
in the spatial structure of plant communities. It is proved theoretically that 
generalist symbionts have only a limited effect on the spatial structure of 
their habitat [23]. In contrast, the inclusion of specialist parasites promotes 
coexistence of the hosts in terms of the existence of a stationary distribution 
under which the density of each host type is positive, while the analysis of 
the corresponding mean-field model supported by numerical simulations sug- 
gests that in any dimension the inclusion of specialist mutualists translates 
into a clustering of the host environment [13]. 

Similarly to most spatial epidemic models such as the contact process, 
the state space of the stochastic processes introduced in [13, 21-23] indi- 
cates whether hosts are either healthy or infected, but does not distinguish 
between different levels of infection of the hosts. However, it is known from 
past research that the number of symbiont individuals, including ectosym- 
bionts, that is, symbionts living on their hosts or in their skin, associated 
to a single host individual may vary significantly. Mooring and Samuel [25] 
found for instance an average of 1791 individuals of the species Dermacentor 
albipictus, commonly known as Winter Tick, on individual elk in Alberta, 
while some individual moose have been found with more than 50,000 ticks. 
In addition, symbionts are generally much smaller organisms than their hosts 
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and reproduce much faster and in greater number. This motivates the devel- 
opment of spatially explicit multiscale models of host-symbiont interactions 
that describe the presence of symbionts through a level of infection of the 
hosts rather than binary random variables (infected versus healthy hosts) 
and include both inter-host symbiont dynamics and intra-host symbiont dy- 
namics. 

In diversity ecology, the infrapopulation refers to all the parasites of one 
species in a single individual host, while the metapopulation refers to all the 
parasites of one species in the host population. In systems involving multiple 
species of parasites, all the parasites of all species in a single individual host 
and in an entire host population are called infracommunity and component 
community, respectively. This terminology shall be employed in this article 
for symbionts in general, that is parasites, commensalists and mutualists, 
even though, strictly speaking, it only applies to parasites. Our main objec- 
tive is to deduce from the microscopic evolution rules of the symbionts, de- 
scribed by transmission rates and reproduction rates, the long-term behavior 
of the metapopulation in a single-species invasion model, and the long-term 
behavior of the component community in a multispecies competition model. 
Since a host species and a symbiont species involved in a symbiotic rela- 
tionship usually evolve at very different time scales (symbionts reproduce 
much faster than their hosts), we shall assume in both invasion and compe- 
tition models that the discrete habitat of hosts is static. This habitat will be 
modeled by a realization of the infinite percolation cluster of supercritical 
site percolation [14]. We shall also assume that symbionts can only survive 
when associated with a host (obligate relationship), which restricts their 
habitat to the infinite percolation cluster, and, to understand the role of 
space on the persistence of the symbiotic relationship, that symbionts can 
only transmit to nearby hosts, adding to the complexity of the interactions. 
In the single-species model, infrapopulations will evolve according to the 
logistic growth process, and the entire metapopulation according to a mix- 
ture of this model and its spatial analog, the contact process [17]. In the 
multispecies model, we will assume that infracommunities evolve according 
to the Moran model [26] , and the entire component community according to 
a mixture of this model and its spatial analog, the voter model [8, 19]. Our 
analysis shows that the condition for survival of a metapopulation strongly 
depends on the carrying capacity of each infrapopulation. Exact calculations 
of the critical curve as a function of the reproduction and transmission rates 
are given when infrapopulations can be arbitrarily large which, as mentioned 
above, is a realistic biological assumption in many symbiotic relationships. 
In systems involving multiple symbiont species, the long-term behavior of 
the component community depends on the spatial dimension: the commu- 
nity clusters in two dimensions whereas coexistence is possible in higher 
dimensions. 
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2. Models and results. The models are constructed in two steps. First, 
the static random environment of hosts is fixed from a reahzation of the in- 
finite percolation cluster of site percolation [14]. This random environment 
naturally induces a random graph. The symbionts are then introduced into 
this universe where they evolve according to an interacting particle system 
on the random graph. The interactions are modeled based on two of the sim- 
plest particle systems: the contact process [17] and the voter model [8, 19]. 
The structure of the random graph implies that the infrapopulation dynam- 
ics are described by logistic growth processes, that is, contact processes on 
a complete graph, and the infracommunity dynamics by Moran models, that 
is, voter models on a complete graph. 

Host environment. To define the habitat of hosts, we set p € (0, 1] and 
let cj be a realization of the site percolation process with parameter p on the 
(i-dimensional regular lattice Z*^, that is, each site of the lattice is either per- 
manently occupied by an individual host with probability p or permanently 
empty with probability I — p. Let EI(a;) denote the set of open/occupied 
sites. By convention, elements of Z'^ and processes with state space S <ZTL'^ 
will be denoted in the following by capital Latin letters. We say that there 
is an open path between site X and site Y if there exists a sequence of sites 
X = Xq,Xi, . . . , Xn = Y such that the following two conditions hold: 

1. For i = 0, 1, . . . , n, we have Xi € ]HI(a;), that is, site Xi is open. 

2. For 2 = 0, 1, . . . , n — 1, we have Xi ~ Xj+i , 

where Xi ~ Xi^i means that the Euclidean norm \\Xi — Xi^i\\ = 1. Writing 
X ^Y the event that sites X and Y are connected by an open path, we 
observe that the binary relation ^ is an equivalence relation on the random 
set EI(cj) thus inducing a partition of E[(cij). In dimensions d>2, there exists 
a critical value pc G (0, 1) that depends on d such that if p > pc then M{u}) 
contains a unique infinite open cluster. The infinite open cluster is also 
called infinite percolation cluster and is denoted by Coo(w) later. We assume 
that p> Pc from now on. Sometimes, the infinite percolation cluster will be 
identified with the graph with vertex set Coo(w) obtained by drawing an 
edge between sites of the cluster at Euclidean distance 1 from each other. 
For more details about site percolation, we refer the reader to Grimmett [14]. 

Random graph structure. In order to define the state space and dynamics 
of the stochastic processes, we first define a random graph Ti^uj) as follows. 
Vertices of are to be interpreted as possible locations for the symbionts, 
while edges indicate how symbionts interact. Let A'^ be an integer and Kjv = 
{1,2,..., N}. The vertex set of 7i{u}) is 

CAr(u;) = {{X,i):X G Coo(w) and i £ K^}. 
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By convention, elements of and processes with state space Cn{uj) will be 
denoted by small Latin letters. To define the edge set, we also introduce 

7r:CAr(u;) — > Coo{(jj) defined by 7r(a;) = X for all x = (X, i) G CAr(a;). 

That is, tt{x) is the Coo(w)-coordinate of vertex x. Let x,y £ Ci\f{uj). Then 
vertices x and y are connected by an edge if and only if one of the following 
two cases occurs: 

1. If tt{x) = iT{y), then x and y are connected by a vertical edge: we write 
x^y. It is convenient to assume that each vertex is connected to itself 
by a vertical edge. 

2. If tt{x) ~ TT{y), then x and y are connected by a horizontal edge: we write 
X y. 

In words, a complete graph with N vertices (which are connected to them- 
selves) is attached to each site of the infinite percolation cluster. Edges of 
these complete graphs are said to be vertical while, for any two sites of the 
infinite percolation cluster, vertices of the corresponding complete graphs 
are connected by edges which are said to be horizontal. Vertical and hori- 
zontal edges correspond, respectively, to potential reproduction events and 
transmission events of the symbionts. 

Invasion of a single symbiont — contact process. To understand the con- 
ditions for survival of a single symbiont species, we introduce a generaliza- 
tion of the contact process [17] on the infinite random graph 'H{oj). This 
defines a continuous-time Markov process whose state space consists of the 
set of the spatial configurations rj:CNioj) — > {0,1}, and whose dynamics 
are described by the Markov generator Li defined on the set of the cylinder 
functions by 

Lif{ri)= [/fe.o) - /(r/)] 

+ E (^E^(y) + Ardegvr(x)E^(?^)) 
x[/(r?.,i) -/(??)], 

where deg7r(a;) is the degree of Tr{x) as a site of the cluster Coo(w), and 
where r]x^i is the configuration obtained from r] by assigning the value i to 
vertex x. Note that the degree of each site of the infinite percolation cluster 
is at least 1, therefore the dynamics are well defined. Thinking of vertices 
in state as uninfected and vertices in state 1 as infected by a symbiont, 
the expression of the Markov generator above indicates that symbionts die 
independently of each other at rate 1, reproduce within their host at the 
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(a) (b) 

Fig. 1. (a) Realization of site percolation with parameter p = 0.6 on the 80 x 80 torus, 
with black squares referring to open sites, and white squares to closed sites, (b) Snapshot 
of the invasion model on the percolation structure starting with a single infected host at the 
center of the universe (a = l and /3 = 2). Each site is represented by a 5 x 5 square, that 
is, complete graph with N = 25 vertices. White squares refer to empty sites, that is, sites 
which are not occupied by a host, black dots refer to symbionts, and gray dots to empty 
vertices. 

reproduction rate a, and transmit their offspring to the nearby hosts at 
the transmission rate /3. That is, each symbiont gives birth at rate a to an 
offspring which is then sent to a vertex chosen uniformly at random from 
the parent's host. If the vertex is uninfected, then it becomes infected while 
if it is already infected then the birth is suppressed. Similarly, each symbiont 
gives birth at rate /3 to an offspring which is then sent to a vertex chosen 
uniformly at random from the hosts adjacent to the parent's host, which 
results as previously in an additional infection if and only if the vertex is 
not already infected. See Figure 1 for simulation pictures of this contact 
process. To study the single-species model, we will sometimes consider the 
stochastic process 

%(X)= r]tix) for anXeCoo(w), 

n{x)=X 

where the sum is over the vertices x € Cn{uj) such that 7r(x) = X . That 
is, fit{X) keeps track of the level of infection of the host at X. This de- 
fines a Markov process whose state space consists of the functions that 
map Coo(w) into {0, 1, . . . , A^} and whose dynamics are described by 

Lif{v)= Yl v{X)[f{f]x^)-fm 
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deg(X) 



x[f{vx+)-fm 



where the configurations fjx- and fix+ are obtained from the configuration f] 
by, respectively, removing and adding a symbiont at site X. In view of the 
geometry of the graph ^.{uj), the stochastic process {fjt}t can be seen as 
a mixture of the contact process with infection parameter /3 on the infinite 
percolation cluster and logistic growth processes with parameter a. 

To describe the predictions based on the invasion model, we let 6i be the 
measure that concentrates on the "all i" configuration restricted to Cn{uj), 
that is, 



We denote by ft the upper invariant measure of the process {r]t}t, which is 
also the limit starting from the measure 5i since the process is attractive. 
The process or metapopulation is said to survive whenever fi ^ 6q and is 
said to die out otherwise. 

First, we observe that, starting with a single infection at time 0, the num- 
ber of symbionts in the system is dominated stochastically by the number 
of individuals in a birth and death process with birth parameter b = a + (3 
and death parameter 1. Recurrence of one-dimensional symmetric random 
walks implies that such a process eventually dies out when 6 < 1 . It follows 
that {r]t}t dies out for all values of whenever q + /3 < 1. 

To find a general condition for survival of the infection, we now assume 
that iV = 1 so that the value of the reproduction rate a becomes irrelevant, 
and compare the process with the one-dimensional contact process. Let T 
be an arbitrary infinite self-avoiding path in the infinite percolation clus- 
ter Coo(w). Since for all sites X eT the degree of X ranges from 2 to 2d, the 
process restricted to the infinite path F, that is, symbionts sent outside F 
are instantaneously killed, dominates stochastically the contact process on F 
with infection parameter /3/d. It follows that the process survives whenever 
/3 > d/3c(l) where /3c (1) is the critical value of the one-dimensional contact 
process, since the self-avoiding path is isomorphic to Z. Standard coupling 
arguments also imply that the survival probability of the infection is non- 
decreasing with respect to both the reproduction rate a and the maximum 
number of symbionts per host A^. It follows directly from these monotonic- 
ity properties that, for all values of N and a, survival occurs whenever 



We now look at the long-term behavior of the metapopulation when N 
is large. As previously explained, this assumption is realistic in a number of 
symbiotic relationships, including the interactions between moose and Win- 
ter Ticks [25]. Under this assumption, at least when the number of symbionts 



5i{r]{x) = i} = l 



for all X £ C]\fi^^) = Coo{uj) x K]\f. 



/3>d/3c(l). 
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is not too large, the stochastic process looks locally like a branching ran- 
dom walk on the random graph 'H{uj), namely the process modified so that 
births onto infected vertices are allowed. In the context of large infrapopu- 
lations, global survival of the metapopulation occurs when the reproduction 
rate a > 1 and the transmission rate /3 > 0. This and the comparison with 
a birth and death process imply that, when N is large and the transmis- 
sion rate /3 is small, a situation which is common in parasitic relationships, 
the metapopulation undergoes a phase transition when the reproduction 
rate a approaches 1. Provided the density of the habitat is large enough, 
the phase transition occurs more generally when the sum of the reproduction 
and transmission rates approaches 1. These results are summarized in the 
following theorem where "survival" means strong survival of the stochastic 
process, that is the existence of a stationary distribution under which the 
density of symbionts is positive. 

Theorem 1 (Contact interactions). Assume thatp>pc and /3 > 0. 

1. For all N > 0, the metapopulation dies out if a + (3 <1 while it survives 
tf(3/d>(3,{l). 

2. If a + (3 > 1 and p is close to 1, then the metapopulation survives for N 
large. 

3. If a + /3/d> 1 and p > pc, then the metapopulation survives for N large. 

As previously explained, the first statement of part 1 follows from a com- 
parison with a two-parameter branching random walk, and the second state- 
ment from a comparison with the contact process restricted to a self-avoiding 
path embedded in the infinite percolation cluster. The proof of the second 
part relies on the combination of random walk estimates and block con- 
structions to compare the process view under suitable space and time scales 
with oriented percolation, and we refer to Section 3 for more details. Sur- 
vival when a > 1 and p > pc in the presence of large infrapopulations can be 
proved based on estimates for the extinction time of the logistic growth pro- 
cess and a new block construction. However, the third part, which indicates 
survival under the weaker assumption a + P/d> 1, can be directly deduced 
from the proof of the second part by again looking at the process restricted 
to an infinite self-avoiding path of hosts. Let F be an infinite self-avoiding 
path, which exists almost surely under the assumption p > pc, and observe 
that, since the degree of each site along this path ranges from 2 to 2d, the 
process restricted to F dominates stochastically the one-dimensional process 
with parameters a and (3/d. The latter survives ii a + P/d> 1 since, under 
this assumption, the proof of the second part indicates that survival occurs 
when p = 1 in any dimension, including d=l. The third part of the theorem 
clearly follows. 
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Competition among multiple symbionts — voter model. To study the in- 
teractions among multiple symbiont species, we introduce the analog of the 
previous model replacing contact interactions with voter interactions [8, 19]. 
The state at time t is now ■.Cn{uj) — > {1,2}, that is, each vertex is occu- 
pied by a symbiont of one of two types. Letting for i = 1,2 

fi{x) = card{y : 7r(y) = 7r(x) and ^t{y) = i}/N, 

gi{x) =card{y:7r(y)~7r(x) and ^t{y) = i}/(iVdeg7r(x)) 

denote the fraction of type i symbionts at site 7r{x) and its neighborhood, 
respectively, the evolution is described by the Markov generator L2 defined 
on the set of the cylinder functions by 

L2f{0= y~] — f t \ , — ^{^1 \^ ^a^^\ \ I a r~^[fi^x,i)-f{0] 

Z;^, .aifi{x) + a2f2(x) + Pigi{x) + I32g2{x) 

, sr^ «2/2(x) +/?2g2(a;) )-f(n] 

,"l/l(^) + «2/2(x)+/3m(x)+/3252(^)^-^^^"''^ 

where (,x,i is the configuration obtained from ^ by assigning the value i to 
vertex x. The transition rates indicate that, regardless of its type, each sym- 
biont dies at rate 1 and gets instantaneously replaced by a symbiont whose 
type is chosen from the nearby symbionts according to the relative fecundi- 
ties and transmissibilities of the two symbiont species. In the neutral case 
when the reproduction rates are both equal to say a and the transmission 
rates are both equal to say /?, the local evolution reduces to the following: 
the type of each symbiont is updated at rate 1 and the new type is cho- 
sen uniformly at random from the same host with probability a/(a + /3) or 
a nearby host with probability /3/(a + /?). Note that the process {Ct}t can 
again be seen as a mixture of two well-known processes, namely, the Moran 
model with selection, and its spatial analog, the biased voter model [6, 7] on 
the infinite percolation cluster. See Figure 2 for simulation pictures of this 
voter model. 

To state our results for the competition model, we set 6 G (0, 1) and denote 
by TT0 the product measure restricted to Cn{oj) = Coo{u^) x Ktv defined by 

Trg{^{x) = 1} = e and 7r4^(a;) = 2} = 1 - 6* for ah 2: G Cat (w). 

From now on, we assume that {£,t}t starts from the product measure irg and 
let =^ stand for convergence in distribution. The process is said to cluster if 

there exists a G (0, 1) such that S^t =^ oSi + (1 — 0)62 as t — )• 00. 

In particular, we have 

lim P{^t{x) ^ ^tiy)) = for all x,y £ Cjv(c^). 
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Fig. 2. Snapshots at time 100 (a) and at time 1000 (b), respectively, of the neutral 
competition model on site percolation with parameter p = 0.6 starting from a Bernoulli 
product measure with density 1/2 (ai = 0:2 = 1/2 and /3i = /32 = l/2j. Each site of the 
lattice is represented by a 5 x 5 square, that is, complete graph with N = 25 vertices. 
White squares refer to empty sites, that is, sites which are not occupied by a host, and 
black and gray dots to symbionts of type 1 and 2, respectively. 

The process is said to coexist if in contrast =^ as t — )• 00 for some vg 
such that 

T^e{i{x) /^(y)} for all x,y^ CN{^^),Xy^y. 

Type 1 is said to invade type 2 if 

p( lim iVt = 00 I A^o = 1) > where Nt = card{x G Cn{^^) ■ ^t{x) = 1}, 

indicating that, starting with a single symbiont of type 1 in the infinite perco- 
lation cluster, there is a positive probability that the number of type 1 keeps 
growing indefinitely. Finally, type 1 is said to outcompete type 2 whenever 
we have the stronger condition =^ 5i . 

Theorem 2 (Voter interactions). Assume thatp>pc. // ai = 02 and 
Pi = 132, the component community clusters in two dimensions, whereas co- 
existence occurs in higher dimensions. 

The analysis of the neutral competition model relies on duality techniques. 
We show, in the neutral case, that the process is dual to a certain system 
of coalescing random walks evolving on the random graph induced by the 
infinite percolation cluster. The long-term behavior of the process is related 
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to the so-called finite/infinite collision property of the graph, which is studied 
in details in two dimensions and higher dimensions separately in Section 4. 

Theorem 3 (Voter interactions with selection). Assume that p = l. 

1- If (^1 ^ Oi2 CLnd Pi > /32, then type 1 invades type 2. 
2. If a2 = and /3i > [32, then type 1 outcompetes type 2. 

Note that, in contrast with the neutral case, when one symbiont type 
has a selective advantage, the transition rates of the competition model 
are no longer linear functions of the local frequencies of symbiont types. 
This leads a priori to a failure of duality techniques. The first part of the 
theorem is established by comparing the number of type 1 symbionts with 
a gambler's ruin model whereas the second part relies on the analysis of 
a certain system of branching coalescing random walks which is dual to 
a particle system related to the model with selection when 02 = 0. Before 
going into the proofs, we note that, though the condition 02 = may appear 
biologically unrealistic, the second part obviously holds when A'' = 1 since 
in this case the value of the reproduction parameters becomes unimportant. 
In particular, the result in the second part is relevant for species in which 
only one symbiont individual can associate with a host individual. Finally, 
we point out that the techniques to prove Theorem 3 also apply to the case 
when p > pc but lead to conditions on the parameters which are far from 
being optimal and to very tedious calculations that only make the key ideas 
unclear. Therefore, for simplicity, we focus on the case p = l only. 

3. Proof of Theorem 1. This section is devoted to the analysis of the 
invasion model, and more precisely to the proof of the second part of The- 
orem 1. The key idea is to show that the branching random walk restricted 
to a large square persist an arbitrary long time provided a + /3 > 1 and is 
large. The combination of our estimates with a block construction implies 
survival of the metapopulation restricted to an infinite self-avoiding path 
of large squares fully occupied by hosts. It is also proved that such a path 
exists whenever the parameter p is close enough to 1. 

The first step is to prove branching random walk estimates in order to 
establish the result in any dimension when p = l, that is, the infinite percola- 
tion cluster consists of the entire lattice. To begin with, we observe that, for 
all M > and 5 £ (0, 1), the process {f]t}t dominates, for sufficiently large, 
the process {Ct}t whose dynamics are described by the Markov generator 

Dif{o= E cw[/(Cx-)-/(c)] 
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+ 



E 



(1 - 5)mX) < M} (aCiX) + A ^ ^^Y)) 



x[/(Cx+)-/(C)], 

where the configurations Cx~ and Cx+ are obtained from ( by, respectively, 
removing and adding a symbiont at site X. Indeed, it suffices that N > M/5 
since in that case 

1 _ jAT-i > 1 _ SjAr^ > (1 - 6)l{j £ [0,M]} for all j = 0, 1, . . . , A^. 

To see this, we observe that the process {Ct}t is a truncated branching 
random walk that allows at most M + 1 particles per site at the same time. 
See Figure 3 where we compare the reproduction rates to site y of a particle 
living at site X for the processes and {Ct}t- on the x-axis we have the 
number of particles at Y, and the parameter x is equal to a if X = y and /? 
if X ~ y. 

Let 6 > such that {a + (3){l — 6) > 1. We will prove, following [4], that for 
all M sufficiently large, the truncated branching random walk {Ct}t survives 
with positive probability, by looking at the branching random walk {Ct}t 
whose dynamics are described by 



Dif{C)= E C{x)[fiCx-)-f{0] 



X [/(Cx+) - /(C)] 

starting with one particle at the origin, where a = (1 — S)a and /3 = (1 
(5)/3/2d. 




M UlS N 

Fig. 3. Reproduction rates for {Ct}t and (* and o, resp.). 
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Lemma 4. For X^O, we haveK[(niX) \ Co(0) = l]>l forn large enough. 

Proof. We observe that E,{Q(X)) satisfies the differential equation 
(see [3], Section 4) 

jMCtix)) = -mm + aKCctix)) + /3 ^ m{Y)), 

whose solution is 

°° on—kj.n 

(1) mix)) =Y,Y1 /^^""'^(O' ^)^^^^"*' 

n=0 A:=0 

where /^("''^^(O, X) is the number of paths from site to site X of length n 
with k loops. To estimate the right-hand side of (1), we let {Uk}k be the 
discrete-time random walk with 

p(TT _ 7 I rr _ - + fory = 2', 

and observe that, for any site X ~ 0, 

y (0, X) " ^ - =P{Un = X\Uo = 0) 

for a suitable Ci = Ci(a,/3) > 0. We refer to [29], Corollary 13.11, for the 
asymptotic estimates of the n-step probabilities. In particular, for X ~ 
and t = n, we obtain 



E(Cn(X)) > j;/i("''=)(0,X) 



n! 
fc=o 

^ ' ^(a + 2d/3)" n! 



^/i("''=)(0,X)- 



V2Tm ^ '(a + 2d/3)" 

> aj^C^n-/^ 

\/27r?i 

for a suitable C2 > 0. Finally, since {a + 2d/3) = {a + f3){l — 5) > 1, we deduce 
that 

E(c„(x),>ii^^M:c.„-/^>i 

provided n is sufficiently large. □ 
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Following the ideas of Lemma 5.3, Remark 5.2 and Theorem 5.1 in [4], 
and using Lemma 4 above in place of [4], Lemma 5.2, one proves that {Ct}t 
survives when M is sufficiently large, and so does, by stochastic domination, 
the metapopulation when N is large and the density p = l. These ideas are 
developed in more details in the following lemma. 



Lemma 5. // (a + /3)(1 — 5) > 1 then the process {Ct}t survives when M 
is sufficiently large. 



Proof. By additivity of {Ct}t, if ^ ~ then the central limit theorem 
implies that 



lim 



PiCni^)>K\CoiO) = K)-l + <^ 



VVar(C„(X) I Co(0) = 



0, 



where the function ^ is the cumulative distribution function of the standard 
normal. Since n is fixed, it follows that, for all e > 0, 

P{CniX) > K for ah X ~ I Co(0) = K) > 1 - e 

for K sufficiently large. Let {Nt}t be the branching process with birth rate 
a + 2d/3 and death rate zero, which represents the total number of particles 
born up to time t. By the same argument as before, there exists C3 > 1 such 
that 

P{Nn < C3K \No = K)>l-e for all K sufficiently large. 

Since, if M > C3K then {Ct}t and {Ct}t coincide (up to time n) on {Nn < 
C3K}, we have 

P{Cn{X)>K\Co{0) = K)>l-2e. 

In order to get 

(2) P{Cn{X) > for ah X ~ I Co(0) = K)>l-2e 

we need to ensure that from time to time n, in no site the process 
{Ct}t on {Nn < C3K} ever exceeds M particles. By geometric arguments 
(see [4], Step 3, for further details), one proves that it suffices to take 
M > 2HqCsK = C4K where Hq is the number of paths of length n in Z"' 
crossing a fixed vertex. To complete the proof, we couple the process {Ct}t 
with a supercritical 1-dependent oriented site percolation process on Z x Z_|_ 
in a way such that the existence of an infinite cluster implies survival for {Ct}t 
relying on the standard rescaling technique introduced in [5]. Let 



Q = {{z, m) G Z X : z + m is even}. 
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and declare site {z,m) G ^ to be good if the host at site zei is infected by 
at least K symbionts at time m x n, where ei denotes the first unit vector. 
Also, let 

Gm = {z £Z:{z, m) is a good site} 

denote the set of good sites at level m. Then, inequality (2) above implies 
that Qra dominates stochastically the set Wm of wet sites at level m in a 1- 
dependent oriented site percolation process on the lattice Q with parameter 
1 — e and with initial condition Wq C Qq (see Durrett [11] for a complete de- 
scription of oriented percolation). The result then follows by choosing e > 
sufficiently small to make the oriented percolation process supercritical. □ 

Since, in the proof of Lemma 4, we consider only the particles of genera- 
tion n, equation (2) holds if, instead of the process {Ct}t^ we deal with the 
process {C"}* obtained by deleting all the particles of generation n' > n. 
In addition, the process {?/"}t, obtained from the metapopulation model 
by assuming that symbionts sent outside [n,?!]*^ are killed, clearly domi- 
nates {crit- 

Lemma 6. Fix n so that Lemma 4 holds. Then, for all e > 0, 
P{Vn{X) >y/N for all X^^\ fil^{0) > Vn) >l-e 
for all N sufficiently large. 

Proof. This follows from (2) (using {f]^}t instead of {C"}*) choosing 
K = \/iV, from stochastic domination when C^^/N /5 < N, and from the 
monotonicity of {r]"}^. □ 

To deduce the second part of Theorem 1 from the previous lemma, we 
use another block construction in order to compare the evolution of the 
metapopulation along an infinite self-avoiding path with oriented percola- 
tion. To apply successfully Lemma 6, any site within distance n of this 
infinite self-avoiding path must be open. The existence of such a path fol- 
lows by choosing the parameter p close enough to 1. First, we fix n so that 
Lemma 4 holds (recall that n only depends on the reproduction rate a, the 
transmission rate /3, and the spatial dimension d). Then, we fix the param- 
eter e > such that 1 — e is greater than the critical value of n-dependent 
oriented percolation. We prove the result when the density p of hosts satisfies 

p>exp((2n + l)-"'logpc), 

where pc is the critical value of site percolation in d dimensions. We tile the 
d-dimensional regular lattice with cubes of edge length 2n + 1 by setting 

Bq = [-n, n]'^ and Bz = {2n + 1)Z + Bq for all Z G Z^. 
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Fig. 4. Picture of the self-avoiding path P. 



Given a realization iv of the site percolation process with parameter p, we 
call a cube Bz open if all the sites X G Bz are occupied by a host, and 
closed otherwise. Our choice of p implies 

P{Bz is open) = p(2n+i)'' > ah Z G Z"^. 

In particular, there exists almost surely an infinite self-avoiding path of open 
cubes, that is, there exists a self-avoiding path {Zi : i G Z} C Z"^ such that 
cube Bzi is open for all i. From this path of open cubes, we construct an 
infinite self-avoiding path of open sites T = {F^ : z G Z} by including all the 
sites belonging to the straight lines connecting the centers of adjacent cubes, 
as shown in Figure 4 where gray squares refer to closed cubes, and white 
squares to open cubes. By construction: 

1. For all z G Z and aU X G F^ + [-n,n]'^, we have X G Coo{^)- 

2. For all z G Z, we have F^ ~ F^+i. 

Site (z, m) G ^ is now said to be good whenever the host at F^ is infected by 
at least \/iV symbionts at time m x n. As previously, we let denote the 
set of good sites at level m. Then Lemma 6 and the fact that the evolution 
rules of the process are homogeneous in time imply that 

(3) P{{z — 1, m + 1) and (z + 1, m + 1) are good | (z, m) is good) > 1 — e 

for sufficiently large A^. Denoting again by Wm the set of wet sites at level 
m in an re-dependent oriented site percolation process with parameter 1 — e 
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the inequality (3) implies that the processes can be constructed on the same 
probability space in such a way that 

P{Wm C On, for all m > I VFo C Go) = 1. 

Since 1 — e is greater than the critical value of oriented percolation, this 
implies as previously that the metapopulation survives, which completes 
the proof of Theorem 1. 

4. Proof of Theorem 2. This section is devoted to the analysis of the 
competition model under neutrality. The process can be constructed graph- 
ically relying on an idea of Harris [16] from a collection of independent 
Poisson processes. In the neutral case, because the transition rates are lin- 
ear functions of the local frequencies, the graphical representation induces 
a natural duality relationship between the spatial model and a system of 
coalescing random walks on Cn{uj), and Theorem 2 follows from certain 
collision properties of symmetric random walks on the infinite percolation 
cluster. 

Duality with coalescing random walks. To define the dual process of the 
competition model under neutrality, we first construct the process graph- 
ically from collections of independent Poisson processes using an idea of 
Harris [16]. Each vertex x £ Ciy{uj) is equipped with a Poisson process with 
parameter 1. Poisson processes attached to different vertices are indepen- 
dent. At the arrival times of the process at x, we toss a coin with success 
probability a /{a + /3) where a is the common reproduction parameter of 
both symbiont types and /3 the common transmission parameter. If there is 
a success, we choose a vertex uniformly at random from the host at site 7r{x) 
and draw an arrow from this vertex to vertex x. If there is a failure, we choose 
a vertex uniformly at random from one of the hosts adjacent to site Tr{x) 
and draw an arrow from this vertex to vertex x. In view of the geometry of 
the graph and the number of vertices per host, this is equivalent to saying 
that: 

— For any pair of vertices x,y £ Cn{uj) with x ^ ?/, we draw an arrow from y 
to x at the arrival times of an independent Poisson process with parameter 
a/{N{a + /3)). 

- For any pair of vertices x,y £ 0^(00) with x -H- y,we draw an arrow from y 
to X at the arrival times of an independent Poisson process with parameter 
/3/(A^deg7r(a;)(a + /3)). 

In any case, an arrow from vertex y to vertex x indicates that the symbiont 
at X dies and gets instantaneously replaced by a symbiont of the same species 
as the symbiont at y. 
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To define the dual process, we say that there is a path from {y,T — s) to 
(x,r), which corresponds to a dual path from (x,r) to {y,T — s), if there 
are sequences of times and vertices 

So = T — s < si < ■ ■ ■ < Sn+i = T and xq = y,xi, . . . ,Xn = x 

such that the following two conditions hold: 

1. for i = l,2,...,n, there is an arrow from Xi^i to Xi at time Sj and 

2. for i = 0, 1, . . . , n, there is no arrow that points at the segments {xi} x 

The dual process starting at {x,T) is the process defined by 

^s(x, T) = {y £ Cj\f{i^) ■ there is a dual path from (x, T) to {y, T — s)}. 

The dual process starting from a finite set of vertices B C Cn^oj) can be 
defined as well. In this case, the dual process starting at {B,T) is the set- 
valued process defined by 

^siB,T) = {y £ Cn{u) : there is a dual path 

from {x,T) to {y,T — s) for some x G B} 

= {y £ C^ii^) ■ y G Csix, T) for some x G B}. 

The dual process is naturally defined only for < s < T. However, it is con- 
venient to assume that the Poisson processes in the graphical representation 
are defined for negative times so that the dual process can be defined for 
all s > 0. Note that, in view of the graphical representation of the competi- 
tion model in the neutral case, the dual process starting at {x,T) performs 
a continuous-time random walk on the random graph Cn{uj) that makes 
transitions 

j z, for z ^ y at rate Q/(A^(a -|- /3)), 
^ \z, for z o y at rate /?/(A^deg7r(y)(a -I- /5)). 

The dual process starting from a finite set B C Cn{uj) consists of a system of 
card(i?) such random walks, one random walk starting from each vertex in 
the set B. Any two of these random walks evolve independently of each other 
until they intersect when they coalesce. This induces a duality relationship 
between the model and coalescing random walks. We refer the reader to the 
left-hand side of Figure 5 for an example of realization of the dual process 
in the neutral case. 

The reason for introducing the dual process is that it allows us to deduce 
the configuration of the system at the current time based on the configura- 
tion at earlier times, but also how vertices at the current time are correlated, 
by keeping track of the ancestry of each symbiont. In particular, the long- 
term behavior of the competition model (clustering versus coexistence) can 
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1 2 3 4 5 6 1 2 3 4 5 6 

Fig. 5. Dual process in the neutral case and branching random walk (^s{x,T). In both 
pictures, N = 2 and for simplicity we have set p=l and d=l. Time goes up, and arrows 
within the same host are drawn in continuous line, while arrows connecting two adjacent 
hosts are drawn in dashed lines. 

be expressed in terms of collision properties of random walks on the infinite 
percolation cluster through the duality relationship between the model and 
coalescing random walks. We now explain this connection in details, starting 
with some key definitions. 

Let G = {V, E) be an infinite connected graph. We call simple symmetric 
random walk on this graph the continuous-time Markov process {Xt}t with 
state space V that jumps from u to u at rate one if and only if {u,v) E E. 
Note that the embedded Markov chain associated to this Markov process is 
the discrete-time random walk {Af„}„, with transition probabilities 

P{Xn+i = V I Xn = u) = — — - if and only if (n, v) G E. 

deg{u) 

Since the graph G is connected, the process {Xt}t is irreducible so either 
all the vertices of the graph are recurrent, in which case the graph is said 
to be recurrent, or all the vertices are transient, in which case the graph is 
said to be transient. Let Xt and Yj be two independent random walks on 
the graph G. Using again the fact that the graph is connected and the Kol- 
mogorov zero-one law, the probability that the two random walks intersect 
infinitely often, namely 

P(for all t there exists s such that Xt+g = Vj+s) 

is either equal to or 1 regardless of the initial positions of the random 
walks. The graph G is said to have the infinite collision property if the 
previous probability is equal to 1, and it is said to have the finite collision 
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property if the previous probability is equal to 0. Such properties for the 
infinite percolation cluster Coo(w) translate through the duality relationship 
with coalescing random walks into coexistence/clustering of the competition 
model, as shown in the next lemma. 

Lemma 7. We have the following alternative: 

1. Coo{uj) has the infinite collision property and then the process clusters, 
or 

2. Coo(w) has the finite collision property and then coexistence occurs. 

Proof. Let B C Cn{uj) be finite and let Ql = {x £ Cn{^) -Ctix) = i} 
denote the set of vertices occupied by a type i symbiont at time t. By 
duality, 

(4) P(GtinS = 0)=E(l-0)l«"*(^'*)l, 

where 9 is the initial density of type 1. Since the number of particles in 
is a nonincreasing function of t and has a limit, the bounded con- 
vergence theorem implies that the probability on the left-hand side of (4) 
also has a limit as t — t- oo. It follows that the process converges to a station- 
ary distribution. To understand how different vertices are correlated under 
this stationary distribution, we take two vertices x,y £ CAr(aj), x ^ y, and 
consider the projections 

X, = 7r{Ux,T)) and ^ = 7r(6(y, T)). 

Let r be the hitting time of the dual processes, that is, 

T = mi{s>0:Ux,T) = isiy,T)}. 

Note that the processes Xs and Yg evolve individually according to continuous- 
time random walks run at rate q := (3/{a + /3) on the infinite percolation 
cluster Coo(w). They evolve independently of each other until time r when 
they coalesce. We set to = and define inductively 

Si = inf{s > : Xg = Yg}, 

ti = inf{s >Si:Xs^ Ys} = inf{s > : / Xs^ or Ys / J 

for i>l. Note that, if the dual processes coalesce at time s then 

ti = Sj+i = oo for all i > max{j : Sj < s}. 

Also, let Mi denote the total number of jumps during the interval of time 
(sj,tj) of either of the dual processes starting at vertex x or vertex y. Wri- 
ting Pi for the conditional probability given the event that Sj < oo and using 
that each dual process jumps to one of the adjacent hosts at rate q and within 
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each host at rate 1 — q, we obtain the fohowing probabiUty: 

oo 

Pi{T > ti) = Y,Pi{T > U I M,=j)Pi{Mi=j) 

(5) ^Efi-^l-'d-^)^ 



3=0 ^ 

N(3 



NJ i_(i_g)(i_i/iv) 



a + N/3 

Let J = J{x,y,uj) = sup{j:sj < oo}, and note that, on the event that the 
dual processes starting at x and y do not coalesce, J is equal in distribution 
to the number of intersections of two independent random walks on the infi- 
nite percolation cluster starting at 7r(x) and vr(y). In particular, if Coo(w) has 
the infinite collision property and if It denotes the number of intersections 
up to time t of two independent random walks starting at 7r(x) and vr(y) 
then (5) implies 

hm P{CT{x)^CT{y)) < lim P{iT{x,T)^iT{y,T)) 

oo j 

T^-oo 



^OO ^ — ^ -LJ. 

j=0 i=0 

lim yf^^yp(/^ = j) = 



T->-oo 

j=0 

by the bounded convergence theorem since P{It = j) ^ as T — t- oo for 
all j G N. This shows that the process clusters. Alternatively, if the infinite 
percolation cluster has the finite collision property, then J is almost surely 
finite so (5) implies that 

hm P{^t{x) ^ ^T{y)) 

T— 5>oo 

= e{l-9)lim P{iT{x,T)^iT{y,T)) 

T— !>oo 
oo 

= 6(1 -e)Y, lim Piirix, T) / iriy, T)\J = j)P{J = j) 

j=0 

which shows that coexistence occurs. This completes the proof. □ 

Note that the previous lemma easily extends to any connected graph in 
which the degree of each vertex is uniformly bounded. That is, given such 
a graph G = {V,E), the competition model can be naturally defined on the 
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graph Gn with vertex set V x IK^v constructed from G in the same way 
as the graph Gn{u}) is constructed from the infinite percolation cluster. 
Then, the proof of the previous lemma implies that, for all N finite, the 
resulting process clusters when G has the infinite collision property, but 
coexists when G has the finite collision property. 

Transience of the percolation cluster and coexistence. Motivated by Lem- 
ma 7, we now prove that the infinite percolation cluster Coo(w) in dimensions 
d>3 has the finite collision property. This follows from the fact that the 
infinite cluster is transient, a result due to Grimmett, Kesten and Zhang [15], 
and that the degree of each vertex is uniformly bounded. We also answer 
the same questions for the infinite percolation clusters in 2 dimensions since 
the proofs are similar, even though this result will not be used to establish 
clustering of the process. 

Lemma 8. The cluster Coo(w) is recurrent in d=2 and transient in d>3. 

Proof. Both statements follow from the fact that any subgraph of a re- 
current graph is recurrent, and equivalently any supergraph of a transient 
graph is transient. This is proved in details in Doyle and Snell [10] based 
on the analogy between random walks and electrical networks so we only 
give an outline of their proof. The idea is to turn the graph under consid- 
eration into an electrical network in which each edge has unit resistance. 
Then, simple random walks on this graph are recurrent if and only if the 
effective resistance of the resulting electrical network between a given point 
and the points at infinity is infinite, as shown in [10]. In other respects, 
Rayleigh's monotonicity law states that if the resistances of a circuit are 
increased, respectively, decreased, then the effective resistance between any 
two points can only increase, respectively, decrease. In particular, removing 
edges induces an increase of the effective resistance between any two points, 
therefore any subgraph of a recurrent graph is recurrent. Again, we refer 
to [10] for the details. 

With the previous result in hands, the first statement follows directly 
from the fact that the infinite percolation cluster in two dimensions is a sub- 
graph of the two-dimensional lattice, which is recurrent by Polya's theorem. 
Transience of the infinite percolation cluster in higher dimensions has been 
established by Grimmett, Kesten and Zhang [15]. Their proof relies on the 
construction of a transient tree-like graph that can be embedded in the in- 
finite percolation cluster. Note that their result applies to bond percolation 
but relies on geometric properties that are known for site percolation as well, 
so their proof easily extends to our context. □ 

Lemma 9. Let {Xt}t and {Yt}t be two independent random walks run 
at rate 1 on Coo(w) both starting at vertex A, and denote by I{X,Y) the 
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number of their intersections. Then, 

EI{X,Y) = oo ind = 2 and EI{X,Y)<oo ind>3. 

Proof. Since the total rate of jump of both random walks equals 2, 

(/•oo \ />oo 

l{Xt = Yt} dt) = 2 P{Xt = Yt) dt 

/•oo 

(6) =2/ P{Xt = B)P{Yt = B)dt 

/•oo 

= 2 Yl {pM,B)fdt, 

where pt{A, B) = P{Xt = B \ Xq = A). Now, we observe that the probability 
that a random walk follows a given directed path from vertex A to vertex B 
is equal to 1 divided by the product of the degrees of the vertices of this 
path excluding the final vertex B. Similarly, the probability that a random 
walk follows the reverse path from vertex B to vertex A is 1 divided by 
the product of the degrees of the vertices excluding the final vertex A, from 
which we deduce that 

(7) {2d)-'ptiB,A)<ptiA,B)<2dpt{B,A) for all 5 G Coo (t^) 

since 1 < deg(^), deg(i?) < 2d. Therefore, when Coo(w) is recurrent, (6) 
and (7) imply that 

/oo 
Y Pt{A,B)ptiB,A)dt 

/•oo 

= d~^ I P2t{A,A)dt = oo, 
Jo 

whereas when Coo(w) is transient, (6) and (7) imply that 



poo 

EI{X,Y)<M Y Pt{A,B)pt{B,A)dt 

/•oo 

= Ad p2t{A,A)dt < oo. 
Jo 



10 

The result then follows from Lemma 8. □ 



Lemma 9 indicates that P{I{X,Y) < oo) = 1 in dimensions d>3, that 
is, Coo(w) has the finite collision property, which, together with Lemma 7, 
implies that coexistence occurs. However, that the expected number of inter- 
sections is infinite does not imply that the number of intersections is infinite 
with positive probability (with probability 1 by the Kolmogorov zero-one 
law). In fact, it is known that recurrent graphs, even with bounded degree. 
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do not necessarily have the infinite cohision property. This has been proved 
by Krishnapur and Peres [20], looking at the comb lattice, that is the sub- 
graph of 7? obtained by deleting all the horizontal edges off the x-axis. 

Infinite collision property of the percolation cluster. We now prove that 
the infinite percolation cluster Coo(w) has the infinite collision property 
in d = 2, which, by Lemma 7, is equivalent to clustering of the neutral 
competition model in two dimensions. We use the same notation as be- 
fore and let {Xt}t and {Yt}t be two independent continuous-time random 
walks run at rate 1 on the infinite percolation cluster. Let Wt = {Xt,Yt) and 
yVn = {Xn,yn) denote the discrete-time Markov chain on Coci^) x C'oo(w) 
with transition probabilities 



where qn{A,B) denotes the n-step transition probability of the lazy sym- 
metric random walk on the infinite percolation cluster. That is, at each time 
step, one of the two coordinates of Wn is chosen at random with probability 
1/2. This coordinate then moves according to the uniform distribution on 
the neighbors or stands still, both with probability 1/2, while the other co- 
ordinate does not change. Note that, at each step, with probability 1/2, the 
process Wn does not move at all. Note also that the processes and 
n can be coupled in such a way that the sequences of states visited 
by both processes are equal. In particular, invoking in addition the Markov 
property and the Borel-Cantelli lemma, to prove the infinite collision prop- 
erty, it suffices to prove that 



The first key to proving the infinite collision property of the cluster is the 
following theorem, which is the analog of Theorem 1 in [1]. We state the 
result in the general d-dimensional case, though we only deal with the two- 
dimensional case in the rest of this section. 

Theorem 10. Let p > pc- Then, there exist a subset il. of the set of 
the realizations with probability one and a collection of random variables 
{5'a}a6Z'* such that the following holds: 

1. We have 5yi(a;) < oo for each to £Q and A G Coo(w). 

2. There are constants ci, 02,03,04 > such that, for all A,B £ Coo{uj), 



P{Wn+l = iA',B')\Wn = {A,B)) 

= A')1{B = B'} + qi{B, B')1{A = A' 



}) 



P{Xn = yn for some n>l\Wo = {A, B)) = 1. 



qn{A,B) > ciu 



d/2 



exp(-C2|^ - -Bp/n) 

whenever \A — B\\/ 5^(u;) < n, 
exp(— C4IA — B\^ /n) 

whenever 5a(<^) < n. 



(8) 



qn{A,B) < csn 



d/2 
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The proof of Theorem 10 follows the lines of the proof of its analog in [1] 
and only differs in two points: first, we consider a discrete-time lazy random 
walk instead of a continuous-time random walk, and second, processes under 
consideration evolve on the infinite percolation cluster of site percolation 
instead of bond percolation. To prove the sub-Gaussian upper estimate, the 
idea is to use a discrete-time version of [24], Theorem 1.1, and the results 
of [9], Sections 5, 6 and 8, while the proof of the sub-Gaussian lower estimate 
follows closely the strategy of [1]. Note that the choice of a lazy random walk 
is motivated by the fact that one cannot expect the lower bound to hold for 
any time n for a standard simple random walk. This is due to the fact that 
it has period 2. In order to avoid unnecessary complications, we prefer to 
deal with an aperiodic random walk. 

In the sequel, to simplify notation, we write sums starting from (or ending 
at) possibly noninteger real numbers, but it is tacitly understood that one 
must consider their integer part. To prove the infinite collision property, we 
define 

F(n) = X;2-" h) qM,X)qn-,{B,X), 

(l-p)n , , 

Fp{n)= 2~-r.\q,{A,X)qn-,{B,X), 

j=pn 

where /9 G (0, 1/2) and ^, B, X G Coo (w) . 

Lemma 11. Fix ^ G (0, 1/2), G Coo(w) and e > 0. ThenF{n)< 

(1 + e)Fp{n) when n is sufficiently large depending on p, A, B, X and e. 

Proof. By the Hoeffding inequality ([18], Theorem 1), we have 




<2exp(-2n(l/2-p)2). 



Taking n such that 

pn > y/n V SAioj) V Sb{u}) 

>\X - A\V\X - B\V SAiuj) V 5b(w), 

we may use the first inequality in (8). Letting $ denote the cumulative 
distribution function of the standard normal, and also applying the central 
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limit theorem, we obtain 



j=pn 
(l-p)n 

' n \ cr / C2n 



(9) 



k=pn 



> E r h7^«-P^ V 



n — J \ n 

2 



(l-p)n 

>(2ci/n)2exp(-C2/((l-p)p)) 2"" 

> C5n-2(2«>((1 - 2p) V^) - 1) 
for some C5 < oo. To conclude, observe that 

, / F(n)-Fp(n) \ ^ , , 

while the previous estimates (9) imply 

F(n)-F,(n) ^ ^.^^ exp(-2n(l/2-p)^) ^ 

n^oo Fp(n) -n^ooC5$((l-2p)Vn)-l/2 

This completes the proof. □ 

With Theorem 10 and Lemma 11 in hand, we are now ready to prove that 
the infinite percolation cluster has the infinite collision property in the sense 
described above, that is, considering continuous-time random walks run at 
a constant rate, say 1. Our proof relies in addition on an argument of Bar- 
low, Peres and Sousi [2] who studied the number of collisions of discrete-time 
random walks moving simultaneously at each time step. In order to under- 
stand the duality properties of the competition model, we need, in contrast, 
to consider a pair of random walks in which only one walk chosen uniformly 
at random can move while the other walk stands still, thus mimicking the 
evolution of a pair of independent continuous-time random walks. 

Theorem 12. Fix a realization uj. Then, for all A,B £ Coo(w), 
P(card{n : = Vn} = 00 \Wo = {A, B)) = I. 

Proof. Let 7 > to be chosen later, and define 

^fc = E E ^(^'^)' 

n=k |X-yl|V|X-B|<v^ 
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where I{X, n) = 1 if there is a coUision at time n at site X , and = otherwise. 
The first step is to find bounds for the first and second moments of Ik when k 
is large. 

Lower bound: E(/fe) > Clog A; for some constant C > which does not 
depend on A, B and for all k > ki{A,B). First, we fix p G (0,1/2) and 
observe that 



^'^/(X, n) = P^'^ {X^ =y^ = X)=Y,2-^h\ q, (A, X)q^., (B, X) 

'l-p)n . . 

2~-^''^jq,iA,X)q^^,iB,X). 



j=pn 

In the previous sum, j and n — j are larger than pn. Hence, for n > A;, if 

\X - A\V\X - B\<^ and k > p~\SA{i^) V Sb{uj)), 
then j A (n-j) > \X - A\V \X - B\V Sa{^) V Sb{^) so Theorem 10 implies 



on pn 
-C2 ( — H : 

J n- J 



n=k\X-A\V\X-B\<y^ j=pn 
Sxi<^)<7 



X exp 

(l-p)n 



E E^- 

n=k\X-A\\/\X-B\<yrpn j=pn 



SxH<i 

This and the central limit theorem imply that, for k large depending on ^4, i?, 
E^'^(4)>exp(-C2/(l-/9)) 

X ^(2ci/n)2card{X e Coo(u;) : 

(10) 



n=k 



\X - A\V \X - B\ <y/pR, 
Sx{i^) < t}- 

Now, by the ergo die theorem, 

card{X G Cooiuj) :\X - A\ V \X - B\ < 5x(w) < 7} 



lim 

n— >oo 



card{X £ Z"^ :\X - A\V \X - B\ < ^} 
P{XeCoo{io),Sxiio)<-f). 



28 D. BERTACCHI, N. LANCHIER AND F. ZUCCA 

In particular, there exists a constant 5 > that only depends on the perco- 
lation parameter p such that for all 7 and k sufficiently large, we have 

(11) caid{X £ Coo(w) : |X - A| V |X - B| < VP^, Sxito) < 7} > 6pn 
for n>k. By (10) and (11), there exists ki{A,B) large such that 

(12) E^'^(4) > c?exp(-C2/(l - p)) V ^ > C6(logA;2 _ j^g^) ^ ^^^^g^ 

^-^ n 

n=k 

for a suitable > not depending on A, B, and all k > ki{A,B). 

Upper bound: IE(/|) < C(logA;)^ for some constant C < 00 which does not 
depend on A, B and for all k > k2{A,B). First, we observe that, for / > n, 

E^'^(/(X,n)/(y,0) = P^'''{Xn = yn = X,Xi = yi = Y) 

= E^'^I{X,7i)E^'^{Y,l-n) 

from which it follows that 

E^.^(/l)<2f;x; E E E^'^A^,-) 

n=k l=n \X-A\\/\X-B\<^Y:\X-Y\<VT^ 
Sx(t^)<7 S'y(a;)<7 

xE^'^/(y,/-n). 
Since /(Yi, Z — n)I{Y2,l — n) = whenever Yi 7^ we also have 
E E^'^/(y,/-n)l{5y(w)<7}<l. 

Y: \X-Y\<yA^ 

Therefore, by applying Lemma 11 twice with e = 1, we deduce that there 
exists 7 large such that for all k sufficiently large 

^2E E 2 E 2-"( ^)'Z.(A^)gn-,(S,^) 



X 



^+2 E E E 2-^'-"^ 



Sy(c^)<7 



X 



/ — n 
i 



xg,(x,y)g,_„_,(x,y) 
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Observing that in the sums over j and i above, we have 
j A (n-j) >pn>pk> Sa{i^) Ssii^), 
i A {I — n — i) > p{l — n) > / p = ^ > Sxioj) 
for ah k large depending on A, B, X, Theorem 10 imphes that 

/ fc^ oo 

7 



X 



+ ^card{y:m\/r^<|X-y| 

< (m + l)Vl - n} 



p 

l=n+^/pm=0 



i=p{l—n) 



i J i{l -n-i) 
X exp(— C4m^(i~^ + — n — i)^^)) . 



In particular, there exists k2{A,B) large such that 



X 



(13) 



7 , f 5(m + 2)^ci 



< 



X exp(-4c4(77T,+ 1)^) 
^-^ n \ p ^-^ I — n j 

n=k \ i=n4-7/p / 



i+7/p 

for suitable constants C7,Cs,Cg < oo not depending on A, B and all k > 
k2iA,B). 

Let k{A,B) = ki{A,B) V k2{A,B). By (12) and (13) and the Paley- 
Zygmund inequality, 

P{h>{Ce/2)logk\Wo = {A,B)) 
>P{h>E{h)/2\Wo = {A,B)) 
> (E^-^(4))V4E^'^(/|) > Ci/4Cg = c> 
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for k = k{A, B) and where, as Cg and Cg, the constant c > does not depend 
on the starting points of the random walks. Then, we define a sequence of 
stopping times and sites as foUows: we start at ng = and (^O) -^o) = {A, B), 
and for all j > 1 we define inductively 

Uj =nj^i + k{Aj^i,Bj^i) and (^^,5^) = (A'„^,3^„J. 

We say that there is a success at round j > 1 when 

card{ne [n^.i, n^) : Af„ = 3^„} > {Ce/2)logk{Aj.i, Bj^i) 

and observe that, at each round, the success probability is larger than c > 0. 
In particular, the probability mass function of the number of successes up 
to round j > 1 is stochastically larger than a Binomial random variable with 
parameters j and c > 0, from which it follows that the ultimate number of 
successes, thus the ultimate number of collisions, is almost surely infinite. 



As previously explained, clustering of the neutral competition model in 
two dimensions follows from the combination of Lemma 7 and Theorem 12. 

5. Proof of Theorem 3. This section is devoted to the proof of Theorem 3 
and the analysis of the competition model in the presence of selection. In the 
asymmetric case, the main difficulty arises from the fact that the transition 
rates are no longer linear with respect to the local frequencies of each sym- 
biont type, which leads a priori to a failure of duality techniques. In particu- 
lar, the invadability of type 1 is proved in the next subsection by comparing 
directly the forward evolution of the competition model with a gambler's 
ruin model. In contrast, under the additional assumption 02 = 0, extinc- 
tion of the symbionts of type 2 is established by invoking duality techniques 
which are available for what we shall call a threshold version of the compe- 
tition model. With this duality relationship in hands, the result follows as 
in the previous section from random walk estimates. 

Invasion of type 1. In order to prove the first part of Theorem 3, we first 
let ^t{X) be the number of type 1 symbionts in the host at X E Z'^ and set 



□ 






2dN{ai+/3i) 
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Observe that Pt{^) = (lt{X) in the neutral case qi = a2 and Pi = (32, and 
that 



Pt{X),qt{X) 



if and only if ^j{X) = C_t{Y) = 0, foraUy^X, 

1 if and only if ^t(X) = Ct{Y) = N, for ah y ~ X. 

Note also that, since qtiX) can take at most {N + l){2dN + 1) different 
values, 



q- ■.= mf{qt{X):qt{X)€{0,l)}>0 

and 

q+:=snp{qt{X):qt{X)€{0,l)}<l. 

Denote by Nt the number of type 1 symbionts present in the system at 
time t. If the number of symbionts of type 1 at time is finite, then 



Nt + 1, at rate ^{N - Ct{X))ptiX), 

X 

Nt-1, atrate^etW(l-PiW), 



X 

where the sum is over all X G Z"^ such that PtiN) £ (0, 1). Now, we observe 
that in the neutral case when ai = 02 and f3i = /32, the embedded Markov 
chain associated to {Nt}t is the simple symmetric random walk on ab- 
sorbed at [note that on each edge {x,y) the rates of invasion between x 
and y are symmetric], therefore the two rates above are equal and 



XX X 

^ Y,{N - UX))pt{X) + ^ 5] UX)pt{X) 

X X 

Y.P^{X). 



X 

This implies that for all ai and /3i such that a\ + (3i^ 0, and for all config- 
urations 

(14) j^l^u^) - L mwTK) E'?*(^)- 

Note that the expression of qt{X) depends neither on a2 nor on /32 therefore 
the equation above does not hold in the neutral case only. It is always true 
regardless of the choice of the reproduction and transmission rates. Now, 
assume that ai > Q2 and /3i > /32- We want to show that in this case the 
process {Nt}t has a positive drift. We say that: 
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1. site X gZ'^ is bad at time t when qt{X) G (0, 1) and Ct{Y) = N for all 
y ~ X and 

2. site X G Z'^ is good at time t when qt{X) G (0, 1) and ^t{y) 7^ N for some 

Note that if A'^o is finite then at any time t the sets of good and bad sites 
are both finite. The first ingredient to proving the result is to observe that 
for any site either good or bad 



Pt 



(X) = (2daiCt{X)+/3, ^tiY)) 



(15) 



X (^dNax + 2(iiV/3i - (ai - a2)(iV - itiX)) 
-(/3i-/52) 5](iV-6(n) 

while if we assume in addition that X is a good site, then 



Pt{X) > 



2dai^t{X)+f3iY:Y^^UY) 



(16) 



2dN{ai+l3i) 
> qt{X){l — c)~^ where c: 



1 



/3i-/52 



2dN{ai + /3i) 



2dN{ai + Pi 



The second ingredient is to observe that, by definition of the lower bound q 
and upper bound q~^, we have q~qt{Xi) < q'^qt{X2) for all Xi,X2 so 



(17) 



q-il-c)-\t{Xi) + q-^qtiX2) 

<q~qt{Xi)+q+{l-cr\t{X2), 

{q"{l-cy'+q+){qt{Xi)+qt{X2)) 

< {q- + q'^){qt{Xi) + (1 - c)-\t{X2)). 



The third ingredient is to observe that if X is bad then all sites Y ^ X are 
good, so the number of bad sites is at most equal to the number of good 
sites. In particular, letting 



and 



B = {X:Xls bad} 
G = {X:X\s good}. 
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there exists a subset G* C G with card(G*) = card(i?). Then, combining 
(14)-(17) gives 

X bad X good 

> 5]g,(x) + (i-cri J] g,(x) + (i-cri Y qt{x) 

xgb x&g* xeG\G* 

^ ^ XeBUG* X6G\G* 

In particular, we have 

By comparing the previous inequaUty with the transition rates of {Nt}t and 
applying the gambler's ruin formula, we can conclude that, starting with K 
symbionts of type 1, we have 

/ ^ -I- + \ 
Pf lim A^t=oo') > 1- ( — -) >0. 



\t-s>oo / \q (1 — c) ^ + 

This establishes the first part of Theorem 3. 

Extinction of type 2. We now prove that, under the extra assumption 
^2 = 0, type 1 outcompetes type 2, which is the second part of Theorem 3. 
Letting Gn denote the Cartesian product of the regular lattice Z'^ and IKtv 
the expression of the Markov generator reduces to 

^ ai/i(x) + + /3292(a;) 

^ aih{x) + Pigi{x) + p2g2{x) 

Introduce /3 = /32 and k = /3i — /32 > 0. Let 7 = k • {2dN)~^ and observe that 
gi{x)^0 implies that K,gi{x) > k ■ {N degTT{x))~^ = k ■ {2dN)~^ = 
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from which it fohows that 

aifi{x) + /3igi{x) _ aifi{x) + Pgi{x) + Kgi{x) ^ (x) +7 



aifi{x) + I3igi{x) + ^2g2{x) aifi{x) + P + Kgi{x) /3 + 7 



aifi{x) + (3igi{x) + (32g2{x) aifi{x) + (3 + Kgi{x) /3 + 7' 

In particular, it suffices to prove that type 1 symbionts outcompete type 2 
symbionts when the dynamics are described by the new process {Ct}t with 
Markov generator 



W(C)= E ^^#^i{5i(^)/o}[/(c.,i)-/(C)] 
^-^ p + 7 



+ E (§T^l{5i(^)/0} + l{ffi(rr)=0})[/(C.,2)-/(C)], 

where as previously C,x,i is obtained from by assigning the value % to ver- 
tex X and leaving the state of all the other vertices unchanged. The process 
can be constructed graphically from the three collections of independent 
processes introduced in Table 1 in the following manner. First, we equip 
the set Cat with a total order relation. Then, at time T„(2;), we have the 
following alternative: 

1. In the case when Uniz) <p:= /3(/3 + 7)""*^, we select a neighboring vertex 
uniformly at random, say y, using in an obvious manner the uniform 
random variable Vn{z) and the total order on the set Cat, and then draw 
an arrow from vertex y to vertex z. 

2. In the case when Un{z) > p:= [3(13 + 7)"^, we draw a set of 2dN arrows 
starting from each of the vertices adjacent to the host at site Tr{z) and 
pointing at vertex z. 

We call the events in rules 1 and 2 above jumping event and branching 
event, respectively. The type of vertex z is updated at times Tn{z), n> 1, 
with z becoming of type 1 if at least one of the arrows that point at vertex z 
originates from a type 1, and of type 2 otherwise. In particular, if all the 

Table 1 
Graphical representation 

Notation Description Interpretation 

Tn{x) Poisson process with parameter 1 Times of an update at vertex x 

Unix) Uniform random variable on (0, 1) Determining tlie new symbiont type 

Vn(x) Uniform random variable on {1, 2, . . . , 2dN} Determining the new symbiont type 
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neighbors are of type 2 then the new type is 2 while if at least one neighbor 
is of type 1 then the new type is chosen uniformly at random from the 
neighbors with probability p or type 1 with probability 1—p, which produces 
the suitable transition rates. We say that there exists a path from space-time 
point {y,T — s) to point (x,T) if there are sequences of times and vertices 

So = T — s < si < • • • < = T and XQ = y, xi,...,Xn = x, 

such that the following two conditions hold: 

1. for i = 1,2, . . . ,n, there is an arrow from Xi^i to Xj at time Si and 

2. for i = 0,1, . . . ,n, there is no arrow that points at the segments {xi} x 

{Si, Sj+l). 

We define a set-valued process {Cs{x,T)}s by setting for all < s < T 

(18) Csix, T) = {y £ Cn ■ there is a path from (y, T — s) to {x,T)}. 

Note that the process {(^six,T)}s consists of a system of branching coalescing 
random walks in which particles independently jump at rate p and branch 
at rate 1 — p. We refer to the right-hand side of Figure 5 on page 19 for 
a picture. The introduction of the process (18) is motivated by the following 
lemma, which is somewhat reminiscent of the duality relationship between 
the biased voter model and branching coalescing random walks. 

Lemma 13. Assume that (^q{z) = 1 for some z £ (^t{-c,T). Then, C,t{x) = 1- 

Proof. Let z E C,t{x,T) with C,q{z) = 1. Then there is a unique path 
from (-2,0) to (x,T). Using the same notation as in the definition of a path, 
we introduce the jump process 

Xt = Xi for all Si <t < Sj+i and Xt = x. 

From the construction of the process {Ct}t7 we have CtiXt) = l for all 0<t<T. 
The lemma follows immediately by applying the equation at time t = T since 
Xt = x. □ 

The rest of the proof relies on standard random walk estimates supple- 
mented with a rescaling argument similar to the one described in Section 3. 
In short, introducing the spatial regions 

B{X, K) = {z€Cn: tt{z) G X + {-K, K)'^} for X G Z'^ and K G Z+, 

the next objective is to show that, for all e > and T = K^, there exists K 
large such that 

P{Ct{x) = 1 for ah x G B{X, 3K) \ Co(x) = 1 for all x G B{X, K)) 

(19) 

>l-e. 
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In view of Lemma 13, inequality (19) follows directly from the following 
result. 

Lemma 14. Let T = . Then, there exists Cio < oo and 710 > such 
that 

P{CT{x,T)r\B{X,K) = 0) < (:7ioexp(-7ioif) for all x G B{X,3K). 

Proof. Let x E B{X,3K). The idea is to define a random walk {VFs}s 
embedded in the system of branching coalescing random walks and connect- 
ing, with probability close to one, the space-time point {x,T) to a vertex 
in the ball B{X, K) at time 0. The process starts at Wq = x. To define the 
dynamics, we also introduce the projection Wg = 7r(Ws). Then, the random 
walk jumps at each time s such that T — s = TniWs-) when we have the 
following alternative: 

1. Jumping event: when UniWs) < there is an arrow from a vertex, say z, 
to Ws-.- Then, the random walk jumps to vertex z, that is, we set Wg = z. 

2. Branching event: when Un{Ws) > p, there are 2dN arrows that point 
at Ws-. Then, the random walk jumps to one of the tails, chosen ran- 
domly and uniformly, that make the random walk's projection Ws closer 
to the center X of the ball. 

Note that Ws £ Cs{x,T) for all s € (0,T), as desired. Introduce 

Y: = \TTi{Ws) -Tri{X)\ for i = 1, 2, . . . , d, 

where vrj is the projection on the ith axis in Z*^. Since each vertex has at 
most dN neighbors which are closer to the center of the target ball for a total 
of 2dN neighbors, we have 

jimh-' . P{Y:^, = y: + l\ V: > O) < p/{2d) =■. r, 

Ihnh-' . P{Y:^f^ = Y:-1\ Y: > O) > p/{2d) + (1 - p)/d =: I. 
ft— >o 

Therefore, {^g^js is stochastically smaller than the random walk {Zs}s with 
a reflecting boundary at zero and that otherwise jumps to the right at rate r 
and to the left at rate /. Let r denote the first time the random walk Zs 
hits the boundary 0. Since r <l, standard large deviation estimates for the 
Poisson distribution imply that 

(20) P{t > CuK I Zo < 3K) < Cu expi-^^K) 

for suitable constants. Finally, we introduce the reverse asymmetric random 
walk {Zs}s with state space Z that jumps to the right at rate I and to the 
left at rate r. Letting 

Uk = P{Zs = for some s > \ Zq = k) 
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a first-step analysis gives l{uk+i — Uk) = r{uk — Uk-i) and tlien 

k-l 

(1 - a){uk - uq) = (1 - a)'^J{ui - uq) = (1 - a''){ui - uq), 

j=0 

where a = r ■ It is straiglitforward to deduce that Uk = a^. In particular, 
P{Zt >K\t<T)< P{Zt >K\Zo = 0) 

(21) 

< P{Zs = for some s <T \ Zq > K) < uk = . 

In conclusion, recalling the definition of the processes {Y^}s, using the 
stochastic domination mentioned above, and applying (20) and (21), we 
obtain 

P{CT{x,T)nB{X,K) = 0) 

< P{Wt i B{X, K))<dx P{Y^ > K) 
<dx P{T>T\Zo<3K) + dx P{Zt >K\t<T) 
<dx Ci2 exp(— 7i2i^) + dx 
for all K large. Since a = r ■ l^^ < 1, the lemma follows. □ 

Lemmas 13 and 14 imply that, when viewed under suitable scales, the 
set of space— time boxes which are void of type 2 dominates oriented site 
percolation with parameter 1 — e. This almost produces the second part of 
Theorem 3. The last problem is that oriented site percolation has a posi- 
tive density of unoccupied sites. To prove that there is an in- all-directions 
expanding region which is indeed void of type 2 symbionts, we apply a re- 
sult from Durrett [12] which shows that unoccupied sites do not percolate 
when e is close enough to 0. Since symbionts of either type cannot appear 
spontaneously, once a region is void of one type, this type can only reappear 
in the region through invasion from the outside. This then implies that our 
process has the desired property and completes the proof of the second part 
of Theorem 3. 
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